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Plastic scintillator (PS) embedding wavelength shifting (WLS) fiber is widely used in high energy particle 
physics, as muon taggers, and also in medical physics and other applications. In this work, a simulation package 
is built to evaluate the effects of the diameter and the layout of the optical fiber on the light yield with different 
configurations. The optimal optical configuration was designed based on the simulation and then validated with 
two PS prototypes under certain experimental conditions. In the study, the top veto tracker (TVT) of the JUNO- 
TAO experiment, comprised of 4 layers of 160 strips of PS was designed and evaluated. When a muon tagging 
efficiency of a PS strip is higher than 99%, the threshold is evaluated. The efficiency of 3-layer out of 4-layer of 
TVT will be higher than 99% even with the tagging efficiency of a single strip as low as 97% using a threshold 
of 10 p.e. assuming 40% SiPM PDE. 


Keywords: Plastic scintillator, WLS-fiber, Light yield, Optical transmission performance, Muon tagging efficiency, JUNO- 


TAO. 


I. INTRODUCTION 


The collisions between the primary cosmic rays and the 
earth’s atmosphere will produce a large number of muons[1], 
the average kinetic energy of which at the sea level is several 
GeV[2]. Because of the high energy, large mass, small decel- 
eration and deflection in the electromagnetic field, and small 
bremsstrahlung effect with the atomic nuclear electric field in 
the matter, the muons will have a strong penetration power[3]. 
A sub-system of the muon veto detector with high muon tag- 
ging efficiency is very important to greatly reduce the back- 
ground induced by the cosmic-ray (CR) muons for the exper- 
iments with only limited overburden near the ground, where 
the flux of muon is normally four to seven orders of mag- 
nitude higher than the underground laboratories with large 
overburden, such as Jinping underground laboratory[4], Gran 
Sasso underground laboratory[5] and Canfranc underground 
laboratory[6]. For example, in neutrino experiments[7—9], 
dark matter experiments[10-13], neutrino-less double beta 
decay experiments[14—16], those muon veto systems require 
muon tagging efficiency higher than 99%. At present, the de- 
tectors based on plastic scintillators (PS) have the advantage 
of easy machining[17—21], flexible structure design, efficient 
and stable performance[22, 23]. The PS detectors[24], es- 
pecially with WLS fibers[25—30] and optical photodetectors 
(Multi-anode PMTs or Silicon PhotoMultipliers(SiPMs)), 
were used in OPERA[31], MINOS[32], LHAASO[33, 34] 
and many other experiments[35-—38]. Meanwhile, the PS de- 
tectors have many applications in geological imaging[39—42], 
2 reactor monitoring and other fields[43—45]. 
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The Taishan Antineutrino Observatory (TAO or JUNO- 
TAO) is a satellite experiment of the Jiangmen Underground 
Neutrino Observatory (JUNO)[8, 9]. The main purpose of the 
TAO experiment is to provide a precise neutrino energy ref- 
erence spectrum for JUNO and benchmark measurements for 
the nuclear database. TAO detector system will consist of a 
central detector (CD), an outer shielding and a veto system. 
The CD will be placed at around 30 m from one core of the 
Taishan Nuclear Power Plant. The CD consists of a 2.8-ton 
gadolinium-doped liquid scintillator (LS) filled in a spherical 
acrylic vessel. The gadolinium-doped LS, as a target material, 
reacts with neutrinos from the reactor to measure the neutrino 
energy spectrum[46]. Because TAO only has a limited ver- 
tical overburden of 4 m of concrete, the major backgrounds 
for the TAO experiment are muon spallation products and ac- 
cidental coincidences, mostly due to the natural radioactivity, 
the top veto tracker (TVT) requires to tag muons with an effi- 
ciency of higher than 99%. 


In this paper, a comparison is realized in sec.II between a 
simulation based on Geant4[47-49] and a measurement of a 
prototype of a PS strip with WLS-fiber readout. The light 
yield results of the experiment and simulation are consistent 
for passing through muons. In sec.III, the diameter and the 
layout of the WLS fiber were further checked for higher light 
yield against with the simulation. An optimized design of the 
PS strip with WLS-fiber and SiPM readout is proposed for the 
TVT system of JUNO-TAO with high light yield and muon 
tagging efficiency, which will provide a good reference and 
guidance for the design of PS detector with WLS-fiber. At 
the same time, the reliability of the optimal design is prelim- 
inarily proved by experiments. In sec.IV, with the proposed 
PS strip design, the expected performance of the TAO TVT 
system is demonstrated. Finally, a summary is given in sec.V. 


63 II. PROTOTYPE OF PS STRIP WITH WLS-FIBER AND 
64 SIMULATION 


6s Muons will deposit their energy when they pass through 
es and interact with the surrounding materials, and the process 
67 of muon energy loss is called muon ionization energy loss[3]. 
es The average energy loss per distance (mass thickness) can be 
e described by the Bethe-Bloth formula (1) [3, 50]: 
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71 where K is a constant, z is the unit charge of the incident 
72 muon, Me and c are the electron mass and the speed of light, 
73 respectively. Z and A are the atomic number and mass num- 
74 ber of the passing-through matter. Wmax is the maximum 
7s kinetic energy that can be transferred to an electron when the 
7s muon collides with the atom. J is the average excitation en- 
77 ergy of the matter. 8 is the ratio of the speed of a particle to 
7s the speed of light. y is the lorentz factor. ô is the correction 
79 factor of the density effect of the matter. The above parame- 
so ters are constant for a given matter. 

s1 From the formula (1), the deposited energy of muon in 
the material is related to the energy of muon and the atomic 
ss number of the material. For thin-layer media with an atomic 
4 number less than 20, such as a PS strip, muons almost pass 
ss through in a straight line. Partially lost energy of muon will 
ss be converted into light in a PS strip, and the light will be ex- 
s7 ported by the WLS-fiber, where SiPMs coupled with the fiber 
© » is an effective, convenient, and rapid method to pick up the 
EN] s photons and then convert them into an electrical signal. 
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Fig. 1. Design of a prototype of the PS strip with WLS-fiber (Sub- 
sequently, it is named option 1) 


s A prototype of the PS with WLS-fiber readout is de- 
ə Signed and fabricated as shown in Figure 1[51], named op- 
ve tion 1. Its size is in a dimension of 2m (Length)x0.1 m 
93 (Width) x 0.02 m (Thickness), and four optical fibers with a 
ə diameter of 1mm are used in total. The pink lines repre- 
ə sent the WLS fiber equally spaced, inserted, and filled in the 
æ surface of the PS strip. The length of the straight portion is 
ə 1.9m. The arrangement is symmetrical in both length and 
ə Width directions. Two of the fibers are focused into a single 


ə group which can be coupled with optical sensors. This op- 
100 tion can reduce the number of optical sensors. For example, 
101 four SiPMs (red circular point in the figure) can be used for 
102 each fiber group, or two PMTs can be used for each PS end. 
103 Finally, the PS is wrapped in reflective film (aluminum foil) 
104 except for the pips of the optical fiber to export photons, and 
105 more details can be found in Ref[51]. The current prototype 
tos Was tested with CR muon before design optimization. 


17 ACR muon survey along the length of the PS strip proto- 
tos type with PMT/SiPM is done with the schema shown in Fig- 
109 ure 2 (the arrangement of optical fibers is not shown). The CR 
110 Muon[28, 52] will be selected by the two muon monitors of 
111 Small scintillators on different hitting locations. When the CR 
112 Muon monitors are triggered, the signals of the PS strip will 
3 be recorded too. 9 points equally spaced along the strip were 
114 measured in total relative to the center of the PS. More details 
115 about the experiment can be found in Ref [51]. The PS strip 
116 is fabricated by Beijing Hoton Nuclear Technology Co. Ltd 
117 [53]. The type of fiber is the WLS fiber BCF92 [54, 55]. For 
18 further understanding and optimization, a Monte-Carlo (MC) 
119 Simulation project is set up based on Geant4[47, 49, 52] with 
120 the parameters provided by the manufacturer of the PS, the re- 
121 flective film and WLS-fiber [55, 56]. The simulation project 
122 mainly includes three parts: the first part is the detector ge- 
123 Ometry, and the PS geometry with optical fibers is designed 
124 through a geometric interface. This is why there are vari- 
125 OUS geometric designs in subsequent optimizations. The sec- 
126 ond is the physical process section, which contains a physical 
127 list and the optical processes of optical photons. The phys- 
128 ical list includes such things as ionization, bremsstrahlung, 
122 multiple scattering, pair generation, Compton scattering, and 
130 photoelectric effects. Optical processes include the genera- 
131 tion of Scintillation and Cherenkov light, wavelength shift 
132 effects, Rayleigh scattering, bulk absorption, and boundary 
133 processes. The third is the extraction and analysis of informa- 
134 tion. The PS and PMT/SiPM are set as sensitive areas; The 
135 PS is responsible for obtaining information about the muon, 
136 and the PMT/SiPM is responsible for obtaining information 
137 about the photons hitting. A parameter interface is provided 
138 in the simulation package to set the properties of the mate- 
130 rial, such as the PS attenuation length, the scintillator yield 
140 (refers to the number of photons generated when the energy 
141 deposited in the scintillator is 1 MeV), the reflectivity of re- 
142 flective film. By scanning these parameters, a series of simu- 
143 lation responses, including the light yield (referring to photo- 
144 electron (p.e.) with consideration of the corresponding pho- 
145 ton detection efficiency (PDE)/quantum efficiency (QE) of 
146 SIiPM/PMT) along the PS longitudinal direction, is obtained. 
147 X? analysis is performed with the experimental data. Table 1 
4s presents a list of parameters corresponding to the minimum 
149 ¥?. In the subsequent optimization work, the parameters of 
150 PS, optical fiber and reflective film in table 1 are the same in 
1s2 Simulation. 


iss The comparison of light yield between the measurement 
154 and simulation at different positions is shown in Figure 3(a), 
155 Where the light yield (in p.e.) refers to the average value of 
156 the photoelectron distribution of the selected muons hitting 
157 each position and the X-axis represents the distance from the 
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Fig. 2. Device and principle diagram for CR muon measurement with the PS strip prototype. The CR muon will be selected by the two muon 
monitors of small scintillators on different hitting locations. When the CR muon monitors are triggered, the signals of the PS strip will be 


recorded too. 
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Table 1. The key parameters of the simulation of the PS strip proto- l 
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Fig. 3. (a): Comparison of light yield between measurement and simulation at different muon hitting positions. The horizontal error bar 
is from the dimension of the muon monitors, and the vertical error bar in MC is the statistical error. Due to the different QE of PMTs and 
coupling, light yield at the two ends is not strictly symmetrical. The light yield of SiPM is larger than that of PMT. Based on the advantages of 
compact requirements, robustness, and higher light yield than PMT, SiPM is proposed as a sensor for TVT. (b): When the CR muon monitor 
in the center of PS, a comparison diagram of the experimental energy spectrum and simulation of PMT at one end. 


to the different QE of PMT, the light yield at the two ends is 
not strictly symmetrical. When using SiPM of the same type 
as a sensor, the distribution of the light yield at one end is 
165 Shown here. As can be seen from the figure, the light yield 
of SiPM is larger than that of PMT. Based on the advantages 
of compact requirements, robustness, and higher light yield 
than PMT, SiPM is proposed as a sensor for TVT. In addition 
to the initial comparison shown in Figure 3(a), another simi- 
lar comparison of energy spectra is also obtained between the 


measurement and simulation at center of PS shown in Fig- 


a 


type 
Material Properties Parameters 
Base Polyvinyltoluene 
PS Scintillation yield (photons/MeV) 8000 
Emission peak (nm) 415 
Attenuation length (cm) 200 @ 400 nm 
: Core Polystyrene 
OES Eee Attenuation length (cm) 380 @ 400 nm 
: Base Aluminium 
Reflective filmi Reflectivity 85% 


ure 3(b), except for individual energy points that do not con- 


iss center of the PS strip in the length direction. When using the 
159 PMT as the sensor, the simulation and data are in good agree- 
160 ment within the error. The light yield at both ends should be 
161 symmetrically distributed around the center of the PS. Due 


form to the experiment, others are in good agreement with 
the experiment within the error range. From an experimental 
perspective, each PMT/SiPM has a different QE/PDE, Since 
the following optimization work does not pay attention to the 
177 impact of electronics, in subsequent studies, the QE of each 
172 PMT or the PDE of SiPM will be input with the same value. 


Il. OPTIMIZATION OF PS STRIP LAYOUT 


179 


On the basis of the agreement between experiment and sim- 
ulation, a further study on the configurations of the PS strip 
design is done including the diameter and the arrangement 
of the WLS-fiber for the light yield and tagging efficiency. 
The relationship between the light yield and fiber diameter is 
shown in Figure 4 for 2 PS configurations with the same num- 
186 ber of fibers embedded. Each point refers to the light yield 
of the muon hitting the PS center. The magenta heart repre- 
sents the measurement of the prototype with PMT mentioned 
earlier. The simulated relationship between the light yield 
and the fiber diameter was represented by red lines with a PS 
width of 100 mm and a spacing of approximately 20 mm be- 
tween neighbor fibers. On the other hand, the blue lines repre- 
sented the simulated relationship with a PS width of 200 mm 
and a spacing of around 40 mm between neighbour fibers. 
Most photons are collected through the WLS fiber by the sen- 
sors. Since the area of the optical sensor (Since TAO TVT 
requirement for readout, subsequent optical sensors default 
to SiPMs) is larger than the dimension of the optical fiber, 
there are still some PS scintillation photons directly collected 
200 by the sensors without going through the WLS fiber. The 
two straight lines at the bottom represent the photons directly 
from the PS scintillation to the SiPM, where the light yield 
contributed by the PS itself is basically independent of the 
fiber diameter. 

According to the trend of the plot, the larger the fiber diam- 
eter, or the smaller spacing of the neighbor fibers, the higher 
the light yield. Finally, a 1.5 mm fiber diameter is suggested 
according to the expected response and reasonable cost. 

Following the requirements of the JUNO-TAO TVT sys- 
tem, the PS strip with 20cm width is suggested according to 
the fabrication, electronics, and cost, but there is still more 
than one option proposed following a different strategy as 
shown in Figure 5. The length and thickness of all the op- 
tions are the same as 2m, and 20 mm is a basic requirement. 
The width of option 2, option 3, option 3-1, and option 4 are 
200 mm, while the width of option 1 is only 100mm as the 
measured prototype. The number of fibers is eight to have a 
similar spacing except option 1. The fiber of option 2 is only 
arranged on the front of the PS as doubled of option 1. Option 
3 is similar to option 2, but with four optical fibers installed on 
the back and another four on the front of the PS, respectively. 
The fiber layout of option 3-1 and option 4 are in a different 
arrangement. The fiber of options 3-1 is more uniform than 
22⁄4 that of option 2, but it is not staggered like that of option 4. 
2s For option 1, option 2, and option 3, the length of the straight 
zee portion of the fiber is 1.9m. For option 3-1 and option 4, the 
2237 length of the straight portion of the fiber is 1.5 m. 
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Table 2. Main configurations of different options of PS strip layout 


Configuration option 1 | option 2 | option 3 | option 3-1 | option 4 

PS width (mm) 100 200 200 200 200 

Fiber numbers 4 8 8 8 8 
Fiber diameter (mm) 1 1.5 15 1.5 1.5 
Fiber spacing (mm) 20 20 20 24 22.5 


The PS width and WLS fiber diameter of all the options are 
229 listed in table 2. The performance (mainly on the light yield 
230 and muon tagging efficiency) and differences among the pro- 
231 posals were further evaluated by simulation. It can be divided 
232 into three categories for comparative study. 


228 


233 For Option 2 VS. Option 3: addresses the dependence of 
234 light yield on fber placement. 


235 For Option 2 VS. Option 3-1: addresses the basic depen- 
dence of light yield on the uniformity of the fber ar- 
rangement in the PS. 


236 


237 


233 For Option 2 VS. Option 4: addresses the dependence of 
239 light yield on the fber layout. 


A. Light yield: nPhotons x PDE 


The angular distribution of the CR muon hitting the PS will 
affect the PS response. To check the difference and eliminate 
the additional influence of angular dependence, the muons 
that are vertically and uniformly incident on the PS strip are 
used in the simulation firstly. The results of the number of 
24 photons received by SiPM for option 1-4 are shown in Figure 
6(a). The X-axis is the number of photons (nPhotons, light 
yield = nPhotons x PDE): the sum of the photons collected 
249 by all the SiPMs in a muon event. The figure also shows 
the average number of photons received by all SiPMs for all 
muon events. The black line is the distribution of the pho- 
tons collected by SiPM of option 1. The red and the green 
lines represent the photon number distributions of option 2 
and option 3, which overlap, indicating that the back or front 
location of fibers has no obvious effect on the collection of 
the photons. The recieved photons of the fiber is from the ab- 
sorption and reemission, reflection of the Cerenkov photons 
of PS rather than the directly generated photons. The average 
value of nPhotons in option 2 is almost twice that of option 1. 
Option 4 has more nPhotons than the previous three options, 
which indicates that the arrangement of fibers has a large im- 
22 pact on the light yield. The difference of nPhotons between 
option 4 and option 2 is studied in detail in section IIIB. All 
the plots have the same entries, therefore the maximum height 
of the plots is related to the distribution width. 

Another simulation is done to evaluate the differences be- 
tween option 2, option 3-1, and option 4. The muon generator 
here is updated in energy and angular distribution according 
to the Ref.[57, 58] as a more realistic situation to model the 
270 response of the PS strip. The results are shown in Figure 6(b). 
The black, blue, and red lines represent the total photon num- 
272 ber distribution of option 2, option 3-1, and option 4, respec- 
273 tively. The figure also shows the average value of photons 
274 received by all SiPM of all muon events. Except for the total 
number of photons, the average value of photons of option 4 
is the highest, option 3-1 is the middle, and option 2 is the 
lowest. Option 4 has the best ability to collect photons. 

From the formula (1), it can be known that the deposited 
energy is directly proportional to the track length of muon 
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Fig. 4. Assuming 30% SiPM/PMT PDE/QE; each point refers to the light yield of the muon hitting the PS center. The magenta heart 
represents the measurement of the prototype with PMT mentioned earlier. Two widths of 200 mm and 100 mm for the PS are checked with 
the same length of 2 m and thickness of 0.02 m. The simulated relationship between the light yield and the fiber diameter was represented 
by red lines with a spacing of approximately 20 mm between neighbor fibers. The blue lines represented the simulated relationship with a 
spacing of around 40 mm between neighbour fibers. The same four fibers are used in each case, which is why the light yield of the 200 mm 
wide PS is smaller than that of the 100 mm one. 
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Fig. 5. All evaluated designs of a single PS strip. The number of fibers is eight with a similar spacing except option 1. The fiber of option 2 is 
only arranged on the front of the PS as doubled of option 1. Option 3 is similar to option 2, but with four optical fibers installed on the back 


and another four on the front of the PS, respectively. The fiber layout of option 3-1 and option 4 are in a different arrangement. The fiber of 
options 3-1 is more uniform than that of option 2, but it is not staggered like that of option 4. 


Option 1 (Fiber diameter 1 mm ): Mean = 61.5 2000 
Option 2 (Fiber diameter 1.5 mm): Mean = 123.9 1800 


option2 : Mean = 239 
option3-1 : Mean = 346 
option4 : Mean = 437 


Option 3 (Fiber diameter 1.5 mm): Mean = 124.2 1600 
Option 4 (Fiber diameter 1.5 mm): Mean = 224.4 
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Fig. 6. (a): With vertical incident muon simulation: Distribution of total number of photons received by SiPM for four different options 
of PS strip layout.nPhotons: the number of photons received/collected by all the SiPMs in a muon event, without the consideration of the 
corresponding PDE of SiPM. Light yield = nPhotons x PDE; the average value of nPhotons in option 2 is almost twice that of option 1. 
Option 4 has more nPhotons than the previous three options.(b): With realistic muon simulation: Distribution of the total number of photons 
received by SiPM with different fiber arrangements. The black, blue, and red lines represent the total photon number distribution of option 2, 
option 3-1, and option 4. 
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Fig. 7. Total photon number versus muon track length. The upper left, upper right, and bottom right of the picture represent the two- 
dimensional diagram of total photon number versus muon track length for option 2, option 3-1, and option 4 respectively. Bottom right: 
comparison of the photon number normalized to unit muon track length. 
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passing through the material. Because of the oblique inci- 
dence of muons of the updated muon generator in the second 
simulation, the track length of the muons in the PS strip can 
be either less(represent an edge events) or greater than the 
thickness of the PS strip. The total photon number in Fig- 
ure 6(b) is much higher than that in Figure 6(a) even with the 
same option: the average value of nPhotons has nearly dou- 
bled, while some signals with small amplitude show up. The 
upper left panel of Figure 7 shows the two-dimensional di- 
agram of total photon number versus muon track length for 
option 2, the upper right panel is for option 3-1, and the bot- 
tom left panel is for option 4. The colors of the three panels 
represent the density of events. The total photon number is 
proportional to the track length as expected. When the track 
length is 20 mm (thickness of the PS), the event density is the 
highest due to the maximum flux density when muon is verti- 
cally incident. The bottom right panel of Figure 7 shows the 
photon number per millimeter of muons passing through the 
PS strip of the three options. It is obvious that the average 
value of photons per unit length (mm) of option 4 is 11.23, 
which is nearly twice that of option 2. This concludes that 
option 4 is the most effective one for the light yield response 
to muon when the muon deposits the same amount of energy 
among these options. 


B. Transmission performance 
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Fig. 8. Position distribution diagram of optical scanning of option 2 
and option 4: the two white frames indicate the simulated points of 
11 different positions along the length direction of the PS, respec- 
tively; and the two red frames indicate the simulated points of 10 
different positions along the width direction of the PS, respectively. 
Each point is equidistant. 


Option 2 and option 4 have a similar configuration on 
WLS-fiber number, PS length, width, and thickness, but their 
light yields show an obvious difference. The photon trans- 
mission process is further checked for a better understanding 
of any specific reasons for the influence of fiber arrangement 


310 


31 


312 
313 
31 


P 


315 
316 


31 


N 


318 
319 
320 


32 


322 
323 
324 
325 
326 


32 


e] 


32 


© 


329 
330 


33 


332 
333 
334 
335 
336 
337 
338 
339 
340 


34 


342 
343 
344 
34 


a 


346 
347 
348 
349 
350 


35 


352 
35: 


Q 


354 


355 


356 


357 
35: 


© 


359 


36l 


[e] 


36 


362 


363 


364 


365 


366 


367 


on the light yield for possible further optimization. 

Muon deposits its energy along its track in the PS and emits 
photons in four 7 directions. Meanwhile, the photons in the 
PS will propagate through attenuation, absorption, and re- 
emission effects, so it is not easy to accurately characterize 
the light transmission performance between different options. 
A specified optical survey is done to check the photon refec- 
tion times before absorption by the WLS-fiber of each photon 
generated by PS, where the survey locations are shown in Fig- 
ure 8 trying to cover the center and edge of the PS strip. The 
orange points are the specified locations to generate optical 
photons. At each point, 15000 photons are generated in four 
m directions to mimic the random photons excited by a muon. 
It is obvious that before the photons enter the fiber, more re- 
flection times, more difficult to reach SiPM (inverse relation- 
ship). At the same time, if the number of photons entering 
the fiber is more, the number of photons arriving at SiPM is 
more (proportional relationship). In order to model the com- 
bination of the two factors, an R-value proposed as the total 
number of photons entering the fiber divided by the average 
number of reflections of photons before entering the fiber. 


Figure 9 shows the distribution of R-value of an optical 
survey at different positions under option 2 and option 4. The 
abscissa of each point is the position corresponding to the 
dash # symbol in the figure, the specific location can be seen 
in combination with Figure 8. The ordinate of each point is 
the R-value stands for optical transmission performance. The 
red line is the R-value of option 4, and the black represent 
option 2. The upper left panel of the figure shows the dis- 
tribution of the R-value along the length direction of the PS 
strip when the position is the center of the width of the PS. 
The R-value of option 4 is higher than option 2 in the whole 
length direction. The lower left panel shows the distribution 
of the R-value along the length direction of the PS strip when 
the position in the edge of the width of the PS strip. The R- 
values of the two are nearly the same, except for the edge in 
the length direction, where the R-value of option 2 is greater 
than option 4. The upper right panel reveals the distribution 
of the R-value along the width direction of the PS when the 
position is the center of the length of the PS, where option 
4 and option 2 have a trend of trade-offs. The lower right 
panel reveals the distribution of the R-value along the width 
direction of the PS when the position is the edge of the length 
direction of the PS, the R-value of option 4 is larger than that 
of option 2 at the positions within 30 mm from the middle in 
the width direction of the PS. 

The results show that when the photon is generated at a 
location where the fiber is sparse, the number of photons en- 
tering the optical fiber will decrease sharply, while the times 
of reflections will also increase. Meanwhile, when photons 
is generated in dense places of fibers, the number of photons 
entering the fiber will increase dramatically, while the times 
of reflections will also decrease dramatically. 

Looking at the four figures as a global picture, with op- 
tion 2 configuration, the maximum vertical coordinate of the 
black point is less than 200, but with option 4 configuration, 
the ordinate of the red point can be greater than 200. The R- 
value of option 4 and option 2 has a trend of trade-offs and 
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Fig. 9. The distribution of the ratio of the photon numbers entering the fiber to the average times of reflections at different positions of the 
optical survey. The value of the X-axis of each point is the position corresponding to the dash # symbol in the graph, and the value of the 
Y-axis of each point in the graph represents the R-value. The upper left panel of the figure shows the distribution of the R-value along the 
length direction of the PS strip when the position is the center of the width of the PS. The R-value of option 4 is higher than option 2 in the 
whole length direction. The lower left panel shows the distribution of the R-value along the length direction of the PS strip when the position 
in the edge of the width of the PS strip. The R-values of the two are nearly the same, except for the edge in the length direction, where the 
R-value of option 2 is greater than option 4. The upper right panel reveals the distribution of the R-value along the width direction of the PS 
when the position is the center of the length of the PS, where option 4 and option 2 have a trend of trade-offs. The lower right panel reveals 
the distribution of the R-value along the width direction of the PS when the position is the edge of the length direction of the PS. The results 
show that when the photon is generated at a location where the fiber is sparse, the number of photons entering the optical fiber will decrease 
sharply, while the times of reflections will also increase. Meanwhile, when photons is generated in dense places of fibers, the number of 
photons entering the fiber will increase dramatically, while the times of reflections will also decrease dramatically. 


the overall difference is small from the lower left and upper 379 to obtain the transmission performance of the whole region of 
right panel of Figure 9. The R-value of option 4 and option 2 ss option 4. The same procedure is followed for the black line 
has a big difference from the upper left and lower right panel ss: as for the red line, and the product of the two integral areas 
of Figure 9. The distribution of R values in the PS region has 2 is used as the transmission performance of the entire region 
been get. To obtain a quantitative relationship between the sss of option 2. Finally, the product of the integral area of the 
optical transmission performance of option 2 and option 4 in 3s4 red line is divided by the product of the black line, and it is 
the entire region, an explanation of mathematical integration sss found that the ratio is 2.21. This shows that the transmission 
(R-value) were provided. The integral area under the red and sss performance of the whole area of option 4 is 2.21 times better 
black lines in the upper left and lower right panels of Figure 9 387 than that of option 2. Therefore, it can be understood why the 
is calculated. Then, the integral area of the red line in the up- sss overall photon number of option 4 is more than that of option 
per left panel is multiplied by the area in the lower right panel ss 2. At the same time, the curve shows that the R-value distri- 
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bution range of option 4 is larger than that of option 2, which 
means that the optical export uniformity of option 4 is worse 
than that of option 2. This is the result of uneven optical fiber 
layout. 

In short, it was found that if the optical fiber in the PS mid- 
dle region is in higher density, the final effective light yield 
will be higher. Therefore, to achieve a higher light yield, it is 
possible to consider placing as many optical fibers as possible 
in the middle of the PS. 

So far, in terms of light yield, many simulations and made 
massive comparisons have done, the best configuration has 
been obtained that is option 4 under the current design con- 
ditions. At the same time, to some extent, the reason why it 
is the best option has been quantitatively explained and good 
reference and guidance suggestions for the processing tech- 
nology and structure design were provided. For a PS detec- 
tor, the most important factor, in addition to light yield, is the 
muon tagging efficiency. In the next section, the muon tag- 
ging efficiency and inefficiency are discussed. 


C. Muon tagging efficiency and inefficiency 


To simplify the simulation process, assuming that the PDE 
of SiPM is equal to 30%. The sum of four SiPM outputs at 
both ends can be acquired under different options to further 
study the muon tagging efficiency. Figure 10 shows the two- 
dimensional distribution of light yield at both ends of option 
2, option 3-1, and option 4 in simulation, respectively. The 
abscissa is the sum of the SiPM outputs at one PS end, and 
the ordinate is the sum at the other end. The range of light 
yield at both ends of option 4 is the strongest and most di- 
vergent. If the sums at both ends of an event are greater than 
their respective threshold at the same time, it is considered 
that this muon has been detected. The event count of muon 
hitting the PS strip is noted as Nau, with the event count of 
detected muon noted as N;iq,. The muon tagging efficiency 
is defined by the ratio Nyag/Nay. Here, the environmental 
background is not considered, mainly to eliminate environ- 
mental background interference. Figure 11 shows the rela- 
tionship between the muon tagging efficiency and the thresh- 
old of option 2, option 3-1, and option 4, respectively. The 
black, blue and red represent the relationship between effi- 
ciency and threshold when the PDE of SiPM is 20%, 30% 
and 40%, respectively. As seen from the three figures, when 
PDE is 20% and the threshold is set to 10 p.e., the tagging 
efficiency of option 4 is still higher than 90%, while option 
2 and option 3-1 are both less than 90%. With the increase 
of the threshold, the efficiency of option 2 decreases fastest, 
followed by option 3-1, and option 4 is the slowest. Accord- 
ing to Figure 11, the required threshold can be found when a 
certain tagging efficiency needs to be achieved. 

Figure 12 shows the corresponding threshold when the 
muon tagging efficiency reaches 99% under different PDEs 
of SiPM. When the PDE of SiPM is 40%, for option 2, the 
threshold cannot exceed 1.5 p.e.; for option 3, it cannot ex- 
ceed 2.4 p.e.; however, for option 4, the threshold can be set 
to 3p.e., which greatly reduces the contribution from dark 
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noise of SiPM. So far, option 4 is superior to other options in 
terms of light yield and muon tagging efficiency. Therefore, 
Option 4 is taken as the current optimal configuration. When 
the PDE of SiPM is 40%, The threshold and corresponding 
tagging efficiency of a single PS strip were summarized in 
table 3. 


Table 3. Threshold and corresponding tagging efficiency of single 
PS strip with 40% SiPM PDE. 


Threshold(p.e.) 
Tagging efficiency 


10 
97% 


3.1 
99% 


6.3 
98% 


15 
96% 


19 
95% 


To further verify the simulation, another prototype with op- 
tion 4 was built, the parameters of PS, optical fiber and re- 
flective film are the same with before prototype in simulation. 
Only the PDE of the backend SiPMs is different. 2 SiPM were 
used as sensors and using the oscilloscope to collect data. The 
SiPM is from the K-series MicroK-40035-TSV[59]. The data 
is saved by the lecroy-HDO4104A oscilloscope. Preliminary 
experimental results were obtained. 

Figure 13 shows the experimental measurement and simu- 
lation results of the PS prototype of option 4 when the muon 
monitor is in the PS center. To avoid the effect of the SiPM 
dark noise and environment background. Signals stronger 
than 3 p.e. was analyzed, and the blue line is the energy spec- 
trum of one SiPM from the measurement. The red line is 
the result of the simulation. Every point has an error bar, the 
abscissa is the light yield, and the ordinate is the log stand 
for the normalized event rate. It can be seen from the figure 
that in the spectrum below 12 p.e., the simulated spectrum al- 
most corresponds to the experiment. It was found that the 
most probable signal amplitude for passing through muons is 
around 8 p.e. from data. So for the final configuration used 
for TAO where there are 4 SiPMs at each end of the PS strip, 
the sum of the most probable signal amplitude will be about 
32 p.e. When the spectrum is in the range of 15 p.e. to 25 p.e., 
within the range of error, the event rate of the simulated spec- 
trum is more than that of the experiment. When the spectrum 
is greater than 25p.e., and the event rate of simulation and 
experiment events is relatively small. In the high energy re- 
gion, the energy spectrum does not show good consistency, 
but it does not affect our optimization work. Several main 
factors lead to differences in the energy spectrum. The first 
is that the coupling between SiPM and PS is not perfect and 
the same for simulation and measurement. The second is that 
the widths of the two small plastic modules for CR monitors 
are not the same, which may lead to a deviation of several 
centimeters in the location where muon hits. Resulting in an 
impact on the energy spectrum. 

In any case, the light yield of option 4 is indeed much 
higher than that of option 1 before optimization. This shows 
that our optimization method is effective. According to the 
above analysis and design, The light yield and the muon tag- 
ging efficiency corresponding to option 2, option 3-1, and op- 
tion 4 have been obtained. When the muon tagging efficiency 
reaches 99%, option 4 has the highest threshold limit. 

To study and discuss the inefficiencies by simulation, the 
track length of all muons in the PS has been counted. The 
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Fig. 10. Assuming 30% SiPM PDE. From left to right: two-dimensional distribution of light yield at both ends of option 2, option 3-1 and 
option 4, respectively. The color represents the density of events. From options 2 to 4, the distribution range of light yield gradually widens. 
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Fig. 11. Relationship between the muon tagging efficiency and the threshold of option 2, option 3-1, option 4, respectively. Where the 
threshold represents the sum of the SiPM of each PS end over threshold individually. As the threshold increases, the muon tagging efficiency 
of option 4 decreases the slowest. 


—k- PDE = 20% s2 Figure 14 shows the ratio of the muon events whose track 
=- PDE =30% so length exceeds certain threshold out of the total events under 
Æ- PDE = 40% soa Option 2, option 3-1, and option 4. The trends are almost the 
sos Same and the difference in ratio was within 1%. The muon 
sos generator, the thickness, length, and width of the PS strip are 
so7 the same for the three configurations, and the declining trend 
sos Of the event ratio is basically consistent under the three con- 
so figurations. The black horizontal dotted line shows the pro- 
portion of 99%, from which it can be obtained that the track 
51 length threshold of the corresponding black vertical dotted 
—— —— —— — 512 line is 3mm. There are 1% muon events whose track length 

prong een opion sis in the PS strip is less than 3 mm. Since the thickness of the PS 
si4 is 20mm, it is certain that the 1% muon events are incident 
sis at a large zenith angle and pass through the edge of the PS, 


When the PDE of SiPM is 40%, for option 4, the threshold can be " which is often referred to as the edge event. Then, if building 


set to 3 p.e., which greatly reduces the contribution from dark noise °" a module of PS to stagger the PS up, down, left and right, 
of SiPM. sis these edge events should also leave longer tracks on the next 


s519 layer or the PS next to them, so as to be triggered. Thus, for a 
s52 PS strip, the edge event may not be tagged(detected), but for 
s21 the whole veto module system, it can be tagged. Therefore, 
s2 under the same threshold, the muon tagging efficiency of a 
ses Single layer module is higher than that of a single PS strip. 
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Fig. 12. When the muon tagging efficiency requirement is 99%, the 
corresponding threshold limits for options 2, 3-1, and 4 respectively. 
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Fig. 13. Comparison between the simulation and experiment in 
terms of the energy spectrum for option 4. Assuming 40% SiPM 
PDE in MC. The PDE of the SiPM sample provided by the man- 
ufacturer is 46% at the wavelength of 420 nm when the working 
voltage is 37.2 V. The data is recorded when the CR muon monitor 
is in the center of PS. 
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Fig. 14. Event ratio which the track length of incident muon in PS 
strip exceeds the threshold in the total events under option 2, option 
3-1, and option 4, respectively. 


Based on this point of view, the JUNO-TAO top veto system 
was designed, called top veto tracker (TVT). 


IV. PERFORMANCE OF JUNO-TAO TVT 


Figure 15(a) is the design schema of the TVT. The left 
panel of Figure 15(a) shows the side view of the whole 4- 
layer PS, the seam of the upper layer of PS corresponds to the 
volume of the next layer of PS, which can eliminate the dead 
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Fig. 15. (a): Design of JUNO-TAO TVT. The left panel is a side 
view. The right panel is a top view.(b): Muon tagging efficiency of 
hitting multi-layer PS with that of single. 


space. At the same time, it is considered that the track direc- 
tion of a muon can be determined by this system. Therefore, 
a 4-layer PS is designed, and the gap between neighboring 
layers is 2 cm. (It refers to the gap between the SiPM end of 
one PS and the end of another PS). The right panel of Figure 
15(a) shows the top view of the whole 4-layer PS. Different 
colors represent different layers of PS. Since the central de- 
tector of TAO is a ball, under the condition of saving funds 
and consumables. Our best effort is to cover the ball’s projec- 
tion surface with the area of PS, resulting in the formation of 
a circle by the four layers. 


Figure 15(b) shows the distribution of the muon tagging ef- 
ficiency of hitting multi-layer PS versus that of a single PS. 
The different colors represent the efficiencies of hitting the 
different number of layers of PS. The black line represents 
the muon tagging efficiency of any one out of four layers. 
The red represents any two out of four layers. Because the 
gap between each layer is very small, the red and black lines 
almost coincide a PS strip efficiency. The blue represents the 
muon that hits the whole four layers and is also tagged by the 
four layers of PS. Since each layer has a certain dead space, 
as the efficiency of the PS strip decreases, the tagging effi- 
ciency reduction on all four layers will be most noticeable. 
The green represents any three out of four layers, from which 
this conclusion can be obtained: even if the efficiency of a PS 
strip is as low as 97%, which is equivalent to the threshold of 
10p.e. can be found in Table 3, the muon tagging efficiency 
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of any three out of four layers will still be higher than 99%. 
In addition, by simulating the environment background (in- 
cluding 788U, ?32Th, 4°K chain), the background event rate 
that can pass the 3 p.e. threshold in any two out of four layers 
of PS is 1280 + 40 Hz and the rate that can pass the 3 p.e. 
threshold in any three out of four layers of PS is < 10 Hz. 
Therefore, the condition can be selected wherein any three 
out of the four layers are triggered to tag muon and reject the 
background. 


V. SUMMARY 


There are numerous PS with WLS-fiber options with differ- 
ent configurations. In this paper, through simulation package, 
the differences between different configurations are checked 
in detail, for example, for the influence of optical fiber diam- 
eter, layout and other factors on transmission performance. 
Based on light yield and muon tagging efficiency, the opti- 
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